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This manuscript is a draft of work in progress, meant 
for network distribution only. It will be updated to a for- 
mal preprint when the numerical calculations will be accom- 
plished. In this draft we develop a consistent closure proce- 
dure for the calculation of the scaling exponents f n of the 
nth order correlation functions in fully developed hydrody- 
namic turbulence, starting from first principles. The closure 
procedure is constructed to respect the fundamental rescaling 
symmetry of the Euler equation. The starting point of the 
procedure is an infinite hierarchy of coupled equations that 
are obeyed identically with respect to scaling for any set of 
scaling exponents f n . This hierarchy was discussed in detail in 
a recent public ation [V.S. L'vov a nd I. Procaccia, Phys. Rev. 
E, submitted, diao-dyn/970507015]. The scaling exponents 
in this set of equations cannot be found from power count- 
ing. In this draft we discuss in detail low order non-trivial 
closures of this infinite set of equations, and prove that these 
closures lead to the determination of the scaling exponents 
from solvability conditions. The equations under considera- 
tion after this closure are nonlinear integro-differential equa- 
tions, reflecting the nonlinearity of the original Navier-Stokes 
equations. Nevertheless they have a very special structure 
such that the determination of the scaling exponents requires 
a procedure that is very similar to the solution of linear ho- 
mogeneous equations, in which amplitudes are determined by 
fitting to the boundary conditions in the space of scales. The 
re-normalization scale that is necessary for any anomalous 
scaling appears at this point. The Holder inequalities on the 
scaling exponents select the renormalizaiton scale as the outer 
scale of turbulence L. 



1. INTRODUCTION 

Unquestionably, the calculation of the scaling expo- 
nents that characterize correlation functions in turbu- 
lence is one of the coveted goals of nonlinear statistical 
physics. In a recent paper |jj , hereafter referred to as pa- 
per I, a formal scheme to achieve such a calculation was 
outlined. In the present paper we discuss the first order 
steps of this scheme, identify the mechanism for anoma- 
lous scaling in turbulence, and discuss the series of steps 
available if one wants to improve the numerical values of 
the computed anomalous exponents. The main ideas how 
to achieve the lowest order closure without breaking the 



rescaling symmetry of the Euler equation are formulated 
such that they repeat essentially unchanged in any higher 
order closure; it seems that no new ideas are necessary. 

To set up the developments of this paper, we present 
now a short review of some essential ideas and results. 
Firstly we stress that a dynamical theory of the statis- 
tics of turbulence calls for a convenient transformation of 
variables that removes the effects of kinematic sweeping. 
In our work we use the Bclinicher-L'vov transformation 
H in which the field v(ro, to\r, t) is defined in terms of 
the Eulerian velocity u(r,t): 



where 



v(r , t \r, t) = u[r + p(r , t), t] 



p(ro,t)= / dsu[r Q + p(r ,s),s] 

J t 



(1) 



(2) 



Note that p(r ,i) is precisely the Lagrangian trajectory 
of a fluid particle that is positioned at r at time t = 
to. The field v(ro, ta\r,t) is simply the Eulerian field in 
the frame of reference of a single chosen material point 
p(ro , t) . Next we define the field of velocity differences 
W(r ,to\r,r',i): 

W(r ,to\r,r',t) = v(ro,t \r,t)-v(r ,t \r',t) . (3) 

It was shown that the equation of motion of this field is 
independent of to, and we can omit this label altogether. 

The fundamental statistical quantities in our study 
are the many-time, many-point, "fully- unfused" , n-rank- 
tensor correlation function of the BL velocity differences 
Wj = W(ro\rjr'j,tj). To simplify the notation we 
choose the following short hand notation: 



X; 



{rjM, tj}, Xj = {r^tj}, W, = yV{Xj) , (4) 



T n {r \X x ...X n ) = (Wi 



(5) 



By the term "fully unfused" we mean that all the coordi- 
nates are distinct and all the separations between them 
lie in the inertial range. In particular the 2nd order cor- 
relation function written explicitly is 



^2 ( r o\ri,r 1 \t 1 ;r 2 ,r / 2 ,t 2 ) 

= (W a (r \r i y i ,t 1 )W l3 (r \r 2 ,r' 2 ,t2)) . 



(6) 
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In addition to the n-order correlation functions the statis- 
tical theory calls for the introduction of a similar array 
of response or Green's functions. The most familiar is 
the 2nd order Green's function G a ^(ro\Xi; x 2 ) defined 
by the functional derivative 



G a0 (r o \X 1 ;x 2 ) 



8W a {r \X 1 ) 
Scp^(ro\x 2 ) 



(7) 



In stationary turbulence these quantities depend on t\ — 
t 2 only, and we can denote this time difference as t. 

The simultaneous correlation function T n is obtained 
from J- n when t\ — t 2 ■ ■ ■ — t n . In this limit one 
can use differences of Eulerian velocities, w(r,r',t) = 
u(r',t) — u(r,t) instead of BL- differences, the result is 
the same. In statistically stationary turbulence the equal 
time correlation function is time independent, and we de- 
note it as 



T n (r 1 ,r' 1 ;r 2 ,r' 2 ; . . .;r n ,r' n ) 
(w(r 1 ,r' 1 ,t)w(r 2 ,r' 2 ,t) . . .w{r n ,r' n ,t)} 



(8) 



One expects that when all the separations i?j = |r.j — 
are in the inertial range, r\ <C Ri <C L, the simultaneous 
correlation function is scale invariant in the sense that 



T„(An, Xr[; Xr 2 , Xr 2 ; Xr n , Ar^ 
X Cn T n (n,r' 1 ;r 2 ,r 2 ; . . . ;r n ,r' n ) . 



(9) 



The goal of our work is: the calculation of the exponents 
from first principles. This is first aim of a statistical 
theory of turbulence. 

One could assume that also the time correlation func- 
tions J- n are homogeneous functions of their arguments, 
including the time arguments. It should be stressed that 
this is not the case, and that in the context of turbu- 
lence, when the exponents are anomalous, dynam- 
ical scaling is broken. The time correlation functions 
are "multi-scaling" in the time variables. In Ref. || 
it was shown that the scaling properties of the time 
correlation functions can be parametrized conveniently 
with the help of one scalar function 2(h). To under- 
stand this presentation, consider first the simultaneous 
function T n (ri . . .r' n ). Following the standard ideas of 
multi-fractals ||,^| the simultaneous function can be rep- 
resented as 



T n (r 1 ,r' 1 ...r' n ) = U n J h ™ d»(h) 
xT n ,h(Pi,p'i,...,p») , 



nh+Z(h) 



(10) 



where U is a typical velocity scale, and Greek coordinates 
stand for dimcnsionless (rescaled) coordinates, i.e. 



rj/Rn , p' = r'JR v 



(11) 



In Eq. (|TJJ) we have introduced the "typical scale of sep- 
aration" of the set of coordinates 



(12) 



At this point L is an undetermined renormalization scale. 
At the end of the calculation we will find that it is the 
outer scale of turbulence. The function Z(h) is defined 
as 



Z(h)=3-V(h) 



(13) 



The function Z(h) is related to the scaling exponents 
via the saddle point requirement 



C n — min \nh 

h 



Z(h)} . 



(14) 



This identification stems from the fact that the integral in 
(|l0| ) is computed in the limit R n /L — > via the steepest 
descent method. Neglecting logarithmic corrections one 
finds that T n oc . 

The physical intuition behind the representation ( JlCj ) is 
that there are velocity field configurations that are char- 
acterized by different scaling exponents h. For different 
orders n the main contribution comes from that value of 
h that determines the position of the saddle point in the 
integral ( [To| ) . It is convenient to introduce a dimensional 
quantity T„./j according to 

/ D \ nh+Z(h) 

T n , h (ny 1 ,...y n ) = u n {^fj 

fnfifa, . (15) 

Dimensional quantities of this type will play an impor- 
tant role in our theory. Especially their rescaling prop- 
erties will be used to organize a 2-covariant theory, see 
Section 3. This quantity rescales like 

T n , h (Xr u Xr[, . . .,Xr' n ) = X^+^T^n,^, ...,r' n ). 

(16) 

Below we will use a shorthand notation for such rescaling 
laws, T ?ih -> X nh+z ^T nA . 

The intuition behind the representation ( |l0| ) is ex- 
tended to the time domain. The particular velocity con- 
figurations that are characterized by an exponent h also 
display a typical time scale tn t h which is written as 



R ( L 

tRJl ~u{R 



(17) 



Accordingly we chose a temporal multi-scaling represen- 
tation for the time dependent function 



ft max 

F n (r \X 1 ,...,X n ) = U n J d»(h) 



nh+Z(h) 



X^ r n,/i(^o|-l,-2, ■ ■ -3n) , 



(18) 
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where J- n: h depends on the dimensionless (rescaled) 
space and time coordinates 



(19) 



As before, we introduce the related dimensional quantity 
J- n ,h according to: 



^„(r () |.V, V„i: r" ( ^ 



nh+Z(h) 



xj=" n)J ,(ro|Si,H 2> ...H n ) (20) 

Of course, we require 

that the function ^n,/i(7*o|Ei, . . . S n ) would be identical 
to T n _h{Pi, p'i, ■ • ■ , p^) when its rescaled time arguments 
are all the same. This representation reproduces all the 
scaling laws that were found for time integrations and dif- 
ferentiations ||. We stress that when the multi-scaling 
representation of turbulent structure functions was intro- 
duced for the first time ( j4|] and references therein) it was 
a phenomenological idea, that could be understood as the 
inversion of the data on £„ , Eq. (O) . We will show below 
that in the context of our theory it appears as a result 
of an exact symmetry of the equations of motion. For 
our purposes it turns out easier to compute the function 
Z(h) than to compute the exponents directly. 

In paper I we showed that the nth order correlation 
functions satisfy, in the limit of infinite Reynolds number, 
an exact hierarchy of equation: 



^«M*r, 



, X„ 



(21) 



drj(r u r 2 ,r)F n+1 (r \X,X,X 2 ,...X n ) =0 . 



The vertex function 7 aM<T (ri, r[, r) is defined as 

r^iru ri,f) = \{[P atl {ri - r) - P afi (r[ - r)} J- 
2 ar„ 



(22) 



The projection operator P a/3 was defined explicitly in 
paper I. It should be stressed at this point that the equa- 
tions of motion do not contain a dissipative term since we 
choose to deal with fully unfused correlation functions. 
For such quantities the viscous term becomes negligible 
in the limit of Reynolds number Re— > oo or the viscos- 
ity v — > 0. This is the main advantage of working with 
unfused correlation functions; the more commonly used 
structure functions do not enjoy this property, and in 
their equations of motion the viscous term remain rele- 
vant also in the limit Re— * oo. The price of working with 
unfused correlation functions is that we have to deal with 
functions of many variables. 

The aim of the present paper is to introduce a system- 
atic method of solving this chain of equations. We begin 
in section 2 with the introduction of the central idea that 



this chain of equations can solved on an "/i-slice" . The 
set of equations for the correlation functions has to be 
supplemented by an associated chain of equations for the 
Green's functions, as is explained in Section 2. In Sec- 
tion 3 we introduce the iJ-covariant closures which pro- 
vide a way to search solutions which are automatically 
scale invariant. In Section 4 we explain in what sense 
Z[h) is obtained from a solvability condition. In Section 
5 we discuss what are the missing calculational steps and 
summarize the present state of the work. 



2. FUNDAMENTAL EQUATIONS ON AN 
"J/-SLICE" 

In this section we introduce equations on an "/i-slice" 
and consider their symmetry under rescaling. Firstly we 
discuss the equations for correlation functions, and sec- 
ondly equations for Green's functions. 



A. Correlation functions and the rescaling group 

Examining the equations (|2l]) one realizes that they 
are linear in the set of correlation functions !F n with 

n = 2, 3, In light of the representation ( p"8| ) of T n in 

terms of J- n ,h we can rewrite Eqs. ( fjlj ) in the form 

jf W dn(h){^-^ n , h (r \X 1 ,...,X n ) (23) 

dr-f(r l ,r 2 ,f)T n+ i,h(ro\X,X,X 2 ,...X n )^ = . 

Since we integrate over a positive measure, the equations 
are satisfied only if the terms in curly parentheses vanish. 
In other words, we will seek solutions to the equations 

—T n ,h(r \X 1 ,...,X n )+ I dff(n,r 2 ,f) (24) 
xF n+1 , h (r \X,X,X 2 ,...X n ) =0 . 

We refer to these equations as the equations on an "h- 
slice". It is important to analyze their properties under 
rescaling. To this aim recall that the Euler equation is 
invariant to rescaling according to 



At- 



(25) 



for any value of A and h. On the basis of this alone one 
could guess that that Eqs. ( ^4| ) are invariant under the 
rescaling group 

ft ► Ar^ , ti ► A 1 ti , J- n ji ► \ nh ZF nj /j . (26) 

In fact we can see that Eqs. ( p4] ) are invariant to a broader 
rescaling group, i.e. 



Xn 



u -» \ l - h u 



^nh-\-2.(h) ^ 



(27) 
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with an n independent function Z(h). This extra freedom 
is an exact result of the structure of the equations (24). 
It is worthwhile to reiterate that this function appeared 
originally in the phenomenology of turbulence as an ad- 
hoc model of multi-fractal turbulence. At this point it 
appears as a nontrivial and exact property of the chain 
of of equations of the statistical theory of turbulence. 

We will show below that the preservation of this rescal- 
ing symmetry will lead to a theory in which power count- 
ing plays no role. As a consequence the information 
about the scaling exponents £„ is obtainable only from 
the solvability conditions of this equation. In other 
words, the information is buried in coefficients rather 
than in power counting. The spatial derivative in the 
vertex on the RHS brings down the unknown function 
2(h), and its calculation will be an integral part of the 
computation of the exponents. It will serve the role of a 
generalized eigenvalue of the theory. 

Of course, we cannot consider the hierarchy of equa- 
tions ( p4j ) in its entirety. We need to find ways to close 
this equation, and this is the main subject of section 3. 
The main idea in choosing an appropriate closure is to 
preserve the essential rescaling symmetry of the prob- 
lem. We will argue bellow that there are many different 
possible closures, but most of them do not preserve this 
rescaling symmetry. We will introduce the notion of Z 
covariance, and demand that the closure does not destroy 
the /i-slice rescaling symmetry (|27|). 



B. Temporal multi-scaling in the Green's functions 

In addition to the n-order correlation functions the sta- 
tistical theory calls for the introduction of a similar array 
of nonlinear response or Green's functions Q m ,n- These 
represent the response of the direct product of m BL- 
velocity differences to n perturbations. In particular 

g 2) ir ii,i 2 ;i3 = ( 7T7 — i — s 

\ dd>(r \x 3 ) 
g^ 1 (r Q \X 1 ,X 2 ,X z ;x i ) 

'3 [W(ro|Xi)W(r |X 2 )W(r |X 3 )] 
54>(r \x 4 ) 

Note that the Green's function G of Eq. (|?|) is Qi t i in 
this notation. The set of Green's functions Q n ,\ satis- 
fies a hierarchy of equations that in the limit of infinite 
Reynolds number is written as 







g*0-^( rQ \X 1 ,X 2 , ■ ■ .,X n ; x n +i) 



dh 

+ J dry^fri.ri.f) 

x g£f i y^ (roXli X 1 ,X 2 , . . .,X n ;x n+1 ) 

= g ( °l a0 - u (r o \X 1 , X 2 . . . X n ,r n+1 ,h + 0)5(h - t n+1 ) 



(28) 



The bare Green's function of (n, 1) order on the RHS of 
this equation are the following decomposition: 

g { °^-^(r \X 1 ,X 2 ...X n ,r n+1 ,t 1 +0) 

= G ^(r i y i ,r n+1 )^Z^(r \X 2) ... ) X n _ 1 ) . (29) 

Note that these functions serve as the initial conditions 
for Eqs. ( |2§| ) at times t n+ i — t±. 

The form of these equations is very close to the hi- 
erarchy satisfied by the correlation functions, and it is 
advantageous to use a similar temporal multi-scaling rep- 
resentation for the nonlinear Green's functions: 



g n ,i(r Q \X 1 ,X 2 



i AT n , x 



n+l) 



dfi(h) 



xQ n ,i,h(ro\Xi,X 2 , ■ ■ .X n ;x n+ i) 



(30) 



^From the rescaling symmetry of the Euler equa- 
tion we could guess the rescaling properties G n ,i,h 
\( n ~ 1 V l ~ 3 g n l h As before, the equations support a 
broader rescaling symmetry group, 



Ar, t -» X 1 - h t, n ,i, 



x (n-l)h-3+Za(h)g nih 
(31) 



In fact the choice of the scalar function Z G (h) is con- 
strained by the initial conditions on the Green's func- 
tions. ^From Eq. (|29| ) we learn that the Green's functions 
depend on Z(h) via the appearance of the correlation 
functions. This means that Zc(h) and Z(h) are related. 
A simple calculation leads to the conclusions that 



Z G (h) = Z(h) . 



(32) 



In this subsection we displayed explicitly only the hier- 
archy of equations satisfied by Q n ,\- Similar hierarchies 
can be derived for any family of higher order Green's 



function Q m „ with 



2,3.... After introducing the 



multi-scaling representation, we can consider the corre- 
sponding Green's functions on an "/i-slice" , Q m ,n,h and 
show that the equations they satisfy have the rescaling 
symmetry of the Euler equation with a Z(h) freedom. In 
all these equations the initial value terms force the scalar 
function Z(h) to be m-independent. 



3. 2-COVARIANT CLOSURES 

Faced with infinite chains of equations, one needs to 
truncate at a certain order. Simply by truncating one 
obtains a set of equations which is not closed upon itself. 
It is then customary to express the highest order quan- 
tities in the truncated set of equations in terms of lower 
order quantities. This turns the set of equation into a 
nonlinear set. Such an approach is not guaranteed to 
preserve the essential rescaling symmetries (27) and (|3l|) 
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of the equations. In this section we develop a systematic 
method to close the hierarchies of equations for correla- 
tion and Green's functions on an "/i-slice" in a way that 
preserves the rescaling symmetry. 

The lowest order closure involves five equations on an 
"/i-slice" . The first two belong to the J- n hierarchy: 



d_ 

di 



r 

-F 2 ,h(Xi,X 2 ) + J dp*y(n,r 2 ,r) 



(33) 



xF 3<h (X,X,X 2 ) = , 

■^T 3 , h (X u X 2 ,X 3 ) + J dfj(n,r 2 ,f) (34) 

xF ith (X,X,X 2 ,X 3 ) = . 

The next pair of equations belongs to the hierarchy of 
G n ,i- Using the representation j3C| ) in (|2g| ) we derive: 



d f 
Qj-Gi,i,h( X i> x 2) + J dr-f(r 1 ,r2,r) 

*G2Xh{X, X- x 2 ) = G { °\x 1 ;x 2 )5(t 1 - t 2 

d f 

Qj-G2,i,h( x i,X 2 ; x 3 ) + J dr~f(r 1: r 2 ,r) 

xG 3A AX,X,X 2 ;x 3 ) 



(35) 



(36) 



= o. 



Here G h stands for the bare Green's function on an 
"/i-slice" . The fifth needed equation is the first equation 
from the third hierarchy of Green's functions G n ,2,h- 



-Gi, 2 ,h(Xr,x 2 ,x 3 ) + J dr~f(r 1: r 2 ,f) 



_d 

di 



(37) 



These five equations are presented symbolically in 
Fig. H and the symbols are explained in Fig. [j]. We now 
show that in the first step of the closure procedure these 
five equations can be considered as i?(/i)-covariant equa- 
tions for five unknowns. These five objects are the 2nd 
order correlation function J- 2l h, the regular Green's func- 
tion Qi i hi and three types of triple vertices. The vertices 
are introduced in Figs. || and [I| We have in these figures 
three relationships that define the vertices Ah, Bh and 
Ch on an "h-slice" in terms of J- 3> h, G2,i,h an d Gi, 2 ,h- 
Note that there is no notion of perturbation theory here 
- we simply define the three vertices in terms of objects 
that appear in Eqs. ©,© and @. 

Eq. (HH) involves the 4th point correlator J-^h, 
Eq. (||) and Eq. (||) involve G 3 ,i,h and G 2 , 2 ,h- These 
are 4th order objects, and we present them in Figs. || 
terms of all the possible decompositions made of lower 
order objects, and in addition new ("irreducible") con- 
tributions which are defined by these relations. In order 
to have a consistent definition we need to add to this 
game the Green's function Gi. 3 ,h- In the context of the 
4th order objects the irreducible contributions are de- 
noted symbolically as empty squares. There are four of 
them, and we denote them as D 3j i t hi f 2,2,(1, Di, 3 ,h an d 



a. Vertices and correlation functions: 

2 

y = ■< , 

r = v , 




a- 




b. Greens' functions: 



^ = 



,1 — 




£L = 



"0.3 





FIG. 1. The diagrammatic notation of the basic objects of 
the theory. Panel a: the vertex 7 and the correlation functions 
T n with n = 2,3,4. Panel b: the bare Green's function 
(thin line), and the dressed Green's functions Gn,m- Objects 
with only straight tails are identically zero. 
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a. Equations for correlators T a h : 




b. Equations for Greens' functions Q^ h - 




c. Equations for Green's functions Q\^ h 



2 




FIG. 2. Diagrammatic reprezentation of the five equations 
for the lowest nontrivial ,2-covarian closure. 





3'rd order correlator: 






i A h / 


(j = AAA <^ 

\ * 


+ AAA/wV 

\ \ 


Ah 


f B h / 


B h r 


+ -WV\AH 


+ AAA/ ^ 


+ aaaaaaV 


Bh 






+ AAA ^ 


+ AAA V 





FIG. 3. Exact reprezentation of the third-order correlation 
function 



Do,4,fi- The first index stands for the number of wavy 
"tails" and the second index for the number of straight 
tails of the empty square. 

A. Systematic Closure 

The first step of closure consists of discarding the ir- 
reducible empty squares. After doing so, we remain with 
precisely five equations (|33|)- ( p7| ) for five unknown func- 
tions. In the next step of closure we retain the empty 
squares as defined by their relations to the 4th order cor- 
relation and Green's function, and add to the list of equa- 
tions on an "h-slice" the equations of motion for the 4th 
order objects, i.e. T^h, Q3,im, Qi,2,h and Qi j3 ^. In 
total we have at this point nine equations. These equa- 
tions will involve four 5th order objects, i.e. ^Fs,h, G&,x,h, 
Q3,2,h and Qi,z,h- Each of these new objects can be writ- 
ten now in terms of all the contributions that can be made 
from low order objects, plus irreducible 5th order vertices 
that we denote as empty pentagons. The second step of 
closure consists of discarding the empty pentagons. This 
gives us precisely nine equations for nine unknown func- 
tions, i.e. ^2,h, Qi,i,h, A hl Bh, C h , and the four empty 
square vertices Dz,i,h, D 2 ,2,h, Di,3,h and £>o,4,/i- 

The procedure is now clear in its entirety. At the nth 
step of the closure we will discard the n + 3th irreducible 
contributions, and will have precisely the right number 



G 




b. Green's functions Q 1 




FIG. 4. Exact representation of the nonlinear Green's 
functions C/2.1 and <5i,2- 



of equations on an "/i-slice" to solve for the remaining 
unknowns. We should stress that this procedure is not 
perturbative since we solve the exact equations on an 
"/i-slice". Our presentation of nth order objects on the 
"/i-slice" in terms of lower order ones is also exact, it 
just defines at the nth step of the procedure a group 
of (n + 3)th new vertices. These vertices are solved for 
only at the (n + l)th step of the procedure, when the 
(n + 4)th vertices are discarded. It would be only fair to 
say that the idea for this procedure came from a care- 
ful examination of the fully renormalized perturbation 
theory for this problem, [jj]. In that procedure one can 
derive equations for the nth order objects that appear 
symbolically like Figs. [5^,, ||b. etc. In addition, one can 
at each step of the procedure have an infinite expansion 
for the irreducible contributions. Nevertheless, the pro- 
cedure explained above is different; firstly, it is consis- 
tently developed on an "/i-slice" , whereas the renormal- 
ized perturbation theory is done for the standard statis- 
tical objects. Secondly, at no point is there any infinite 
expansion whose convergence properties are dubious. We 
just go through a set of explicit definitions solving an ex- 
act set of equations. The only question that needs to be 
understood is the speed of convergence of this scheme in 
terms of the scalar function Z(h) which parametrizes the 
anomalous behavior. 



B. .2-covariance 

A crucial property of our closure procedure is that it 
guarantees that power counting remains irrelevant on an 
"/i-slice" for an arbitrary step of the procedure, and the 
scalar function Z(K) cannot be computed from power 
counting. To see this we need to find the rescaling prop- 
erties of the triple and higher order vertices. We start 
with the triple vertices Ah, Bh and Ch- The first one is 
defined by its relation to Q1.2J1, see Fig. ^b. Using the 
facts that 

T 2>h - X 2h+Z(h) T 2<h , (38) 
Qi,i,h -» \ z(h) - 3 Gi,i,h , (39) 
Gi,2,h^\- h+z(h) - 6 Gi,2, h ,. (40) 
we find that Ah has to transform according to 

A h - X 2h - 2Z ^- 9 A h . (41) 

Armed with this knowledge we proceed to the definition 
of Bh through its relation to </2,i,ftj equation in Fig. ||a. 
We can check that the rescaling property of 

e 2 , Vl -> X h+Z(h) - 3 Qi,2,h (42) 

agrees exactly with the two terms that contain the vertex 
Ah on the RHS. Accordingly also the term containing Bh 
has to transform in the same way, leading to 

B h -» X 4h - 2Z ( h )- e B h . (43) 
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4.1, 



a. 4 th order correlator: 

j- AAAAAA 1 
i WW J 





b. Green's function £J 3 1 




= 3{ 



AAA 

AAAAAA 



+ 21 




c. Green's function Q 2Xh : 



2,2,/i 




= 3{ 



AAA/ 



+ 9 




FIG. 5. Exact representation of the forth-oredr correlation 
function (Panel a) and nonlinear Green's functions (Panels b 
and c). Empty squares reprezent irreducible contributions 
to the 4th order vertices which are neglected in the lowest 
order iJ-covariant closure. Double line represent either wavy 
or straight line. For more details see Paper I. 



In making this assertion we assumed that there is no can- 
cellation of the leading terms in the equation. Otherwise 
the vertex Bh would be smaller. Lastly we use the def- 
inition of Ch by the relation to see Fig. [|, that 
transforms like 



3,h 



(44) 



All the terms that include Ah and Bh have the same 
rescaling exponent as that of J-^,h- We can therefore 
find the rescaling exponent of Ch'- 



(45) 



We again assumed that there is no cancellation of the 
leading terms in the equation for F^h- If there is can- 
cellation, the vertex Ch can be smaller. 

At this point we can check that the first step in our 
closure scheme leads to a Z{h) covariant procedure. Con- 
sider firstly the three contributions of Gaussian decom- 
position which are first on the RHS of the equation in 
Fig. [|a. These rescale like \ 4h + 2Z ( h ) ( anc l their ratio to 
the LETS is proportional to {R/L) z( - h \ For 2(h) pos- 
itive these contributions become irrelevant in the limit 
(R/L) — > 0. The 45 contributions that come next con- 
tain pairs of triple vertices, and we need to use the rescal- 
ing properties (|4l|), (fl3|), and J45| ) to find their rescaling 
exponents. We find that they all share the same rescal- 
ing exponent . In hindsight this should not be surprising. 
This is a result of the assumption that in the definitions of 
the three vertices there are no cancellations in the leading 
scaling behavior. Thus the rescaling exponent of all the 
45 contributions could be obtained from analyzing one 
of them. The nontrivial fact is that the common rescal- 
ing of all these terms is exactly the rescaling of the LHS 
of the equation, which is \ 4h + z ( h ) . This means that 
our closure for the 4th order correlation functions can- 
not introduce power counting. Note that the rescaling 
neutrality with respect to counting of h and of natural 
numbers follows from the rescaling symmetry of the Eulcr 
equation, and is shared also by Gaussian contributions. 
On the other hand the neutrality with respect to Z{h) 
is nontrivial, and follows from a judicious choice of the 
proposed closure scheme. We will refer this property as 
2-covariance. 

It is important to understand now that the proposed 
closures for the other 4th order objects, like Fig. ||b are 
also iJ-covariant. All that changes is the number of wavy 
and straight tails on the LHS and RHS of the equations, 
and the rescaling exponents change in the same way on 
the two sides of the equation. 

We can now consider the next step of closure, taking 
into account the irreducible 4th order vertices (empty 
squares), discarding the 5th order empty pentagons (ir- 
reducible contributions, for details, see Paper 1 The pro- 
cedure follows verbatim the one described above for the 
triple vertices and the rescaling exponents of the irre- 
ducible 4th order vertices D m n h are determined: 
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D m . nJl - \ d — {h) D m ^ h , (46) 
d m ,n{h) — 2nh — 32(h) — 3m — 4 . 

We can check the terms that appear in the 5th order cor- 
relation and Green's functions, which are made of com- 
binations of triple and 4th order vertices. We discover 
that all these terms share the same rescaling exponent, 
and that it agrees precisely with the rescaling of the 5th 
order correlation and Green's function. Accordingly also 
the second step of closure is Z-covariant. 

It becomes evident that we develop a systematic Z- 
covariant closure scheme, and that power counting will 
not creep in at any step of the procedure. In the next 
section we show that the scalar function Z(h) can be 
computed from this scheme as a generalized eigenvalue. 
We also explain the role of the boundary conditions in 
the space of scales, and how the re- normalization scale is 
chosen. 



condition that the linear operator had an eigenvalue zero, 
and the functions that we seek would have been identified 
as the "zero-modes" associated with this eigenvalue. Our 
equations are not linear, and the solvability condition is 
not that simple. Nevertheless we know apriori that at 
every step of the closure we will need to find the solvabil- 
ity condition of the nonlinear set of equations, and this 
condition will determine the numerical value of Z(h). 

Even after determining Z(h) the statistical functions 
will be determined only up to a factor which may depend 
on h, which is the measure fj,(h) in Eqs([l8|) and (|3"c|). 
In order to determine this factor we will need to fit all 
correlation functions T n to the boundary conditions in 
the space of scales. At that point the computed values of 
Z{h) will determine whether it is the inner or the outer 
scale of turbulence that appears as the renormalizaiton 
scale. We will show in the lowest step of closure that the 
outer scale L is selected. 



THE SCALAR FUNCTION Z AS A 
GENERALIZED EIGENVALUE 



5. SUMMARY AND THE ROAD AHEAD 



In this section we demonstrate that at any step of our 
closure, the scalar function Z(K) can be found only from a 
solvability condition. We will also explain the role of the 
boundary conditions in the space of scales in determining 
the scaling functions in this theory. 

The point is really rather simple. First observe that 
our initial equations for correlation and Green's functions 
on an "/i-slice", like ( p4| ) or (35)-(^) are linear functional 
equations. The equations for the correlation functions 
are not only linear, but also homogeneous. The equa- 
tions for the Green's functions are not all homogeneous, 
but as we explained in the last section the inhomogeneous 
terms are much smaller than the homogeneous terms in 
the limit (R/L) — > 0, and they can be discarded. This 
means, of course, that all the functions on an "/i-slice" 
can be determined only up to an over all numerical con- 
stant. At this point we may have even more than one 
overall free constant; on the face of it the equations for 
the correlation functions are independent of the equa- 
tions for the Green's functions, and we can have different 
overall constants in the correlation and the Green's func- 
tions. 

Every step of closure turns a set of linear functional 
equations into a set of nonlinear functional equations. 
We claim that nevertheless all the functions appearing in 
these equations can be determined only up to an overall 
numerical constant. The extra freedom of many possi- 
ble constants disappears now, since the correlation and 
the Green's functions are coupled after closure. But one 
overall constant remains free. The reason is of course 
the property of .Z-covariance, which includes as a special 
case invariance to an overall scaling factor, or rescaling 
by A . 

If our equations were linear, this freedom would have 
meant that we need to require the standard solvability 



Up to now we described the general ideas and how 
the closure scheme should work. The main point of the 
analysis is that the equations of motion of the statisti- 
cal objects admit an exact rescaling group that contains, 
for a given "/i-slice", one unknown scalar function Z(h) 
whose calculation is sufficient for the evaluation of the 
scaling exponents ( n via Eq. (|l4|). 

Unfortunately, the actual calculations that are called 
for in this scheme are far from trivial. In the lowest 
step of closure we are faced with five coupled integro- 
differential equations in many space-time variables, and 
the complexity increases rapidly in the higher steps of clo- 
sure. We are currently attempting to solve numerically 
the lowest step of closure with the aim of demonstrating 
the existence of anomalous scaling. We are not particu- 
larly interested in the first step in precise values of the 
exponents £„, it is more important to show that they dif- 
fer from their K41 counterparts. Numerical accuracy and 
the convergence of Z(h) will be examined in the higher 
steps of closure. The complexity of the numerics means 
that this is a long program that is expected to last for 
a couple of years. For this reason we chose to describe 
the ideas of the closure scheme even at this stage when 
the numerical implementation is lacking. It is our feel- 
ing that the procedure is sound and that there is a good 
chance to obtain results that will justify the effort that 
is called for in the numerical implementation. 
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